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We show that any quantum information theory based on anticommuting operators must be sup¬ 
plemented by a superselection rule deeply rooted in relativity to establish a reasonable notion of 
entanglement. While quantum information may be encoded in the fermionic Fock space, the unre¬ 
stricted theory has a peculiar feature: The marginals of bipartite pure states need not have identical 
entropies, which leads to an ambiguous definition of entanglement. We solve this problem, by prov¬ 
ing that it is removed by relativity, i.e., by the parity superselection rule that arises from Lorentz 
invariance via the spin-statistics connection. Our results hence unveil a fundamental conceptual 
inseparability of quantum information and the causal structure of relativistic field theory. 
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I. INTRODUCTION 

Over the past two decades, quantum information (QI) 
theory has been developed into a rich and highly success¬ 
ful framework for information processing. Operating in 
the domain of quantum mechanical Hilbert spaces, many 
QI tasks can be approached from a viewpoint that is 
strongly influenced by computer science, while the phys¬ 
ical systems represented by the Hilbert spaces are some¬ 
times of secondary concern. This abstraction from the 
physical context is a significant virtue of QI theory. Deal¬ 
ing with problems purely on the level of a Hilbert space, 
its subsystems, and operations thereon, without refer¬ 
ence to the specihc physical implementation, provides a 
level of freedom and generality that is highly desirable. 
Statements can be made for all Hilbert spaces of a certain 
type. For instance, quantum information processing with 
qubits can typically be investigated without reference to 
their implementation — although examples exist, where 
specifying the encoding of the qubit is relevant for an 
abstract problem, see, e.g., Ref. [1]. Besides multi-qubit 
systems, quantum harmonic oscillators are prominent ex¬ 
amples for successful abstraction. Based on the commu¬ 
tation relations, a bosonic Fock space is constructed, that 
provides a playground for quantum optics, irrespective of 
the particular realization, be it as optical modes, super¬ 
conducting circuits, or vibrational degrees of freedom, to 
name only a few examples. 

Here, another type of Hilbert space — the fermionic 
Fock space — will be considered. That is, the basic alge¬ 
bra is based on anticommuting operators, rather than 
commuting ones. Absent physical interpretation, one 
may yet work with such a Hilbert space, identify its sub¬ 
systems, and their correlations. In other words, one may 
attempt to construct an abstract fermionic QI theory, 
see, e.g., Refs. [2-10]. Conceptually, it is of great im¬ 
portance to collect all types of particles encountered in 
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nature in a common framework, and hence strengthen the 
generality of QI theory as a whole. However, as we shall 
discuss, the physically unrestricted fermionic QI theory 
suffers from a disconcerting malady: As noted already a 
decade ago [11-13], the marginals of bipartite fermionic 
pure states may not have matching spectra. This leaves 
the typical notion of entropy of entanglement in a state 
of ambiguity due to the mismatch of reduced state en¬ 
tropies. Depending on the choice of subsystem, different 
amounts of entanglement would be attributed to the sys¬ 
tem. Indeed, facing a globally pure state, for which one 
subsystem is maximally mixed, while the other remains 
pure would be a significant concern, for instance, in con¬ 
nection with Hawking radiation and the black hole infor¬ 
mation paradox (see, e.g., [14]). These problems do not 
occur in theories with a natural tensor product struc¬ 
ture, like bosonic modes or qubits, where the Schmidt 
decomposition guarantees symmetric marginal entropies 
for pure states. For fermions, on the other hand, map¬ 
pings to a tensor product space, i.e., to qubits, do not 
generally preserve the structure of the subsystems [15], 
and the issue persists. 

In this work, we resolve this problem. We show that it 
is overcome by imposing a superselection rule (SSR) that 
forbids coherent superpositions of even and odd num¬ 
bers of fermions. This property of the SSR, which is 
equivalent to the restriction to parity-conserving opera¬ 
tor (sub-)algebras, is most fortunate, since the SSR also 
gives rise to a natural definition of subsystems [9], and 
is hence widely used. Although, as we shall show, the 
problem of asymmetric pure state marginals is thus re¬ 
moved, it seems rather artificial to enforce such a re¬ 
striction within the abstract theory. In particular, since 
the often referenced argument by Wick, Wightman and 
Wigner [16] in defense of this SSR is based on time re¬ 
versal symmetry, which cannot be an exact symmetry of 
nature in the face of charge-parity (CP) violation [17, 18] 
and the charge-parity-time (CPT) theorem [19-21]. In 
contrast, we discuss the later argument for the SSR 
from [22] invoking invariance under rotations by 27r. Cru- 
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daily, we note that the causal structure of relativistic 
quantum field theory (QFT) enters into this line of rea¬ 
soning, i.e., Lorentz invariance is required to establish 
the spin-statistics connection. We argue that, only once 
the fermionic model is embedded in a physical context, 
in this case relativistic QFT, does the SSR arise natu¬ 
rally. The significance of this approach is twofold. First, 
it ensures that the entropy of entanglement remains a 
well-defined concept in physical QI theories. Second, it 
provides a new perspective on the motivation for using 
the SSR, which has been employed in detailed studies of 
fermionic entanglement, such as [9]. Third, this result 
knits together the fabrics of QI and relativistic QFT in 
a fundamental way: We argue that the fermionic theory, 
and by extension all quantum information theory, must 
be viewed in the physical context of relativity. 

Besides the possible interest for fermionic QI theory 
and applications such as entanglement within QFT in 
curved spacetimes, this work adds a new facet to the dis¬ 
cussion of informational constructions of quantum theory 
(see, e.g.. Refs. [23-25]), by introducing an information- 
theoretic aspect of SSRs — a fascinating topic in its own 
right (see, e.g., Refs. [16, 26-33] for a selection of litera¬ 
ture). It is also of interest to note that the SSR does not 
remove the intrinsically different character of fermionic 
modes and qubits, as indicated by the existence of pure 
states that satisfies the SSR, but still cannot be consis¬ 
tently mapped to multi-qubit states, see the Appendix. 

In the following, we will first outline the construction 
of the fermionic Fock space, as well as of the pure and 
mixed states in such a Hilbert space in Sec. II. To under¬ 
stand the origin of the problem described above, we will 
then discuss the subtleties involved in forming subsys¬ 
tems of fermionic modes in Sec. Ill, and give an example 
for a pure state that features marginals with different 
entropies. Finally, the role and the origin of SSRs are 
discussed in Sec. IV, and we show how the problem can 
be disposed of in Sec. V. 

II. FERMIONIC FOCK SPACE 

Let us consider a system of n fermionic modes with 
mode operators b/, and 6 ^ for (A: = 1,..., n), which satisfy 
the anticommutation relations 

{b,,bj}^0, (1) 

for all i,j. The vacuum state is annihilated by all bk, i.e., 
bk II0)) = 0 V/c, and the purpose of the double-lined no¬ 
tation for the state vectors will become apparent shortly. 
The creation operators bj. populate the vacuum with sin¬ 
gle fermions, that is, ||0)) = || 1^)). When two, or more, 
fermions are created, the corresponding tensor product 
of single-particle states needs to be antisymmetrized due 
to the indistinguishability of the particles. We use the 
convention 

&I^ll|0))= l|U))A||M) = ||U))||M) , (2) 


where we use the double-lined notation to imply the an¬ 
tisymmetrized wedge product “A” between single-mode 
state vectors with particle content (as opposed to the 
standard notation | • ) | • ) = | ■ ) ( 8 ) | • )). With this 
definition at hand, and postponing possible physical re¬ 
strictions, one may write arbitrary pure states on the 
Fock space as 

n 

ll'I^)) = 70 II 0 )) + 5^7. 111.)) + E A.'c 111.)) Illfc)) + ■ • ■ ■ 

2=1 j,k 

(3) 

where the complex coefficients 70 , 7 ^, 'jjk, . ■. are chosen 
such that the state is normalized. Similarly, mixed state 
density operators can be written as convex sums of 
projectors on such pure states. For more details on this 
notation and the fermionic Fock space, see, e.g., [15, 34]. 


III. PARTITIONING THE FERMIONIC FOCK 
SPACE 

The basic ingredients for an abstract fermionic QI the¬ 
ory are pure states on the total Hilbert space, \ip) S 
H, and reduced states with respect to some bipartition 
{A\B), i.e., pA(B) = ATB{A){\'^){tp\). With this, one 
may already study correlation measures such as the mu¬ 
tual information, and, most importantly, the entropy of 
entanglement £{\iIj)) = S{pa) = S{pb), where S{p) = 
—Tr(pln(p)). As the first step in such a construction, 
it is hence necessary to establish a meaningful notion 
of subsystems. Since the particle number in the Fock 
space need not be fixed, we will consider entanglement 
between different modes. However, due to the antisym- 
metrization, the Fock space is not naturally equipped 
with a tensor product structure with respect to the indi¬ 
vidual mode subspaces. These subspaces may nonethe¬ 
less be defined by invoking consistency conditions [15] 
that ensure that the expectation values of all local ob¬ 
servables Oa {i.e., as in, operators pertaining only to the 
modes of the subspace A) yield the same result for the 
global state pab i and for the corresponding local reduced 
states PA = TrB(pyiB), i.e., 

= (^^)p. • (4) 

This procedure uniquely defines the mode subspace 
marginals of any global state, i.e., the partial trace op¬ 
eration in the following way. For an arbitrary number of 
modes, up to a scalar prefactor, all density matrix ele¬ 
ments of the reduced state pA that is obtained by “tracing 
out” one mode labelled by k can be brought to the form 
^pi • ■ • bj,, II 0 )) (( 0 II b,,,... b,,. for some pm,Vn^k where 
m = \,... ,i and n = 1,... ,j. The information about 
these matrix elements arises from ignoring (or “forget¬ 
ting”) whether the mode k is occupied or not. That is. 
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the partial trace operation may be defined as 


T, 4 bl...bl bl\\ 0 )){{ 0 \\bkb.,...b,^) 

= Tvk{bl...bl ||0))((0|| 6.,...6.,) 

= II 0))((0 II61.1 ...61., , ( 5 ) 

while all other elements of pab, proportional to oper¬ 
ators of the form bj^^ .. .bj^. || 0 ))(( 0 || 5^ . .. b^. and 

II 0 )) (( 0 II ... b„- do not contribute to the 
reduced state when k is traced out. The rule can hence be 
expressed as follows: Operators corresponding to modes 
that are being traced out are anticommuted towards the 
vacuum projector before being removed. The ordering of 
the operators on the left-hand side of Eq. (5), that is, the 
position of the operators b^ and b^. relative to the other 
annihilation and creation operators, is determined by the 
condition in Eq. (4). Intuitively, it can be seen that the 
ordering is chosen such that the number of commutations 
to move bk and b\ towards the projector on the vacuum 
state keeps track of any signs incurred when applying ob¬ 
servables Oa to the state pab before the partial trace and 
commuting these past the operators bj, and bj,. A more 
detailed derivation and a rigorous proof of this argument 
is given in Ref. [15]. An equivalent notion of subsystems 
may alternatively be obtained by restricting the operator 
algebras to the corresponding subalgebras for A and B, 
see Ref. [ 10 ]. 

At this point, it is helpful to understand the differ¬ 
ences between fermionic modes and qubits. For any fixed 
number n, the fermionic n-mode Fock space is isomor¬ 
phic to an n-qubit space. A widely known example for 
such an isomorphism is the Jordan-Wigner transforma¬ 
tion, (see, e.g., [9]). Such mappings generally do not 
commute with the procedure of partial tracing [15], since 
local mode operators are mapped to global qubit opera¬ 
tions. In other words, it is generally not possible to estab¬ 
lish isomorphisms between a fermionic n-mode state and 
an n-qubit state in such a way that also all of the respec¬ 
tive fermionic marginals are isomorphic to their qubit 
counterparts. An illustration of this problem is shown 
in Fig. 1. Consequently, the (quantum) correlations be¬ 
tween n fermionic modes may generally not be identified 
with those of the isomorphic n-qubit states. 

In spite of this inequivalence, the partial trace, and 
hence the subsystems and their entropies remain well de¬ 
fined for fermionic modes. Moreover, the construction of 
the density operators and its marginals is based solely 
on the algebraic structure of (1), together with the re¬ 
quirement that the expectation values of subsystem ob¬ 
servables yield the same result when evaluated using ei¬ 
ther the global states or the corresponding marginals. No 
other assumptions are required for a consistent definition 
of the subsystems, and their total correlations. For in¬ 
stance, the mutual information Iab^ a measure of the 
overall correlation between subsystems A and B, is in 
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FIG. 1. Inconsistency between fermions and qubits: 

The global n-mode fermionic state pab may be mapped to 
an isomorphic n-qubit state pab- The marginals oi pab, e.g., 
Pa = Trs(p.4B), are well-defined by Eq. (4), and may also be 
mapped to isomorphic qubit states {e.g., pa Pa)- However, 
as shown in Ref. [15], it is in general impossible to match all 
the marginals of the n-mode state to all marginals of the n- 
qubit state, Pa A PA = TrBipAB). An example for a state 
featuring this problem can be found in the Appendix. 


this context already well defined by the expression 

Iab{pab) = S{pa) + S{pb) — S{pab) , (6) 

where pa(b) = ^t:b(A){pab), and S{p) = -Tr(pln(p)) 
is the von Neumann entropy of the density operator p. 
But, as we shall elaborate on shortly, the same cannot 
be said for genuine quantum correlations, i.e., entan¬ 
glement. Consider, for instance, the non-superselected 
two-mode pure state given by 

11^)) = 70 ||0)) -F7fe ||lfe)) +7fc' l|lfc')) +lkk' ||lfc))||lfc')) • 

(7) 


According to the prescription of Eq. (5), the single-mode 
reduced states can be quickly checked to be 


Pk = Trfe/(]]?/;))((V'll) = (|7o|^ + \lk'?) 

+ (l7fcP + l7fefc'Hl|lfc))((lfcll (8a) 

+ (7o7fc+7fc'7L') l|0))((lfell+H-c. , 

Pfc. =Tr4||r/>))((V^[[) = ([7o|^ + l7fcn l|0))((0|| 

+ (l7fcf+ l7fcfcf)IIM)((U'll (8b) 

+ [(7o7fc' -7fc7fcfe') l|0))((lfe'|l +H. c.l. 


where the symmetry between the subsystems is broken 
by the different relative signs within the off-diagonal el¬ 
ements. The eigenvalues of the two reduced states do 
not match in general. For example, when 70 = 7 fc = 
Ik' = Ikk’ = 1 / 2 , the mode k appears to be in a pure 
state (with eigenvalues 0 and 1 ), whereas the state of 
the mode fc' is maximally mixed (both eigenvalues are 
1/2). Normally, the entropy of the subsystem of a pure 
state would be considered as an entanglement measure. 
Here, depending on the choice of subsystem, one would 
either conclude that the overall state is maximally en¬ 
tangled, or not entangled at all. This problem is not 
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limited to pure states. It persists for mixed states, where 
the entropy of entanglement is of central importance for 
the entanglement of formation. Such an ambiguity in 
the definition of entanglement is of course highly unde¬ 
sirable. One possibility to resolve the issue, would be to 
change the definition of entanglement, and work with a 
non-symmetric quantity. On the other hand, such a dras¬ 
tic step may not be required, if one is willing to embed 
the abstract fermionic quantum information theory in a 
physical framework. As we shall show in the following, a 
reasonable definition of entanglement between fermionic 
modes is obtained when invoking an additional physi¬ 
cal principle—the spin-statistics connection, which itself 
arises from (special) relativity. 


IV. INVOKING 

RELATIVITY—SUPERSELECTION RULES 

No SSRs have been introduced up to this point. Note 
that the term SSR may refer to different restrictions. 
For instance, they may arise from fundamental symme¬ 
tries of the system, such as parity [16], or charge conser¬ 
vation [26-28]. Alternatively, effective (or generalized) 
SSRs originate from practical limitations, such as parti¬ 
cle number conservation due to energy constraints, see, 
e.g., [29-31, 35, 36]. Both type of restrictions may be for¬ 
malized as constraints on the observables, see, e.g., [7, 8 ]. 

Here, we will formulate such constraints in a different, 
but equivalent way, as restrictions on the components 70 , 
7 i, 7 jfe, ... [see Eq. (3)], of pure state decompositions 
with respect to the Fock basis. In particular, we will con¬ 
sider any coherent superpositions of even and odd num¬ 
bers of fermionic excitations as unphysical. The argu¬ 
ment that we will employ to defend this position is based 
on the well-known spin-statistics connection, relating the 
anticommutation relation algebra to the transformation 
properties associated to objects of half-integer spin. Re¬ 
call that, a priori, we have made no assumptions on the 
physical realization of the anticommutation relations of 
Eq. (1). Nonetheless, as an inevitable consequence of this 
anticommutation relation algebra, the excitations of the 
mode operators must satisfy Fermi-Dirac statistics. 

To explain this, we must surrender some level of ab¬ 
straction and provide a physical context. When we place 
the fermionic quantum information theory in the context 
of relativistic field theory, the spin-statistics connection 
follows from Lorentz invariance (see, e.g., Refs. [37, 38] 
and [39, pp. 52[). Thus, the (special) theory of relativity 
enforces that we must interpret the fermionic excitations 
of our theory as particles of half-integer spin. In other 
words, relativity establishes a connection between the 
algebra of the mode operators and the transformation 
properties under spatial rotations of the corresponding 
fermionic particles. In particular, the states of odd num¬ 
bers of these fermions switch their sign under rotations 
by 27r, while even numbers of fermions are left invariant. 
If superpositions of even and odd numbers of fermions 


were permitted, rotations by 27r would change the rela¬ 
tive sign of the respective contributions within the super¬ 
position. For example, suppose that an observer was able 
to prepare a superposition || -|-)) = (|| 0 )) -I- || 1^ ))) /v^ 
in some fermionic mode k. By rotating their frame 
by 27r, the state ||-|-)) is converted to the orthogonal state 
II “)) = (II 0 )) ~ II Ifc )))/v^- Here, one may think of the 
observer physically rotating their preparation device rel¬ 
ative to the detector distinguishing the orthogonal states 
||-|-)) and ||—)). In other words, if one demands that rota¬ 
tions by 27r should leave physics unchanged, one thereby 
enforces the parity SSR that forbids superpositions of 
even and odd numbers of fermions [22]. 

Crucially, this line of reasoning is intimately tied to 
relativity, whereas the argument for the parity SSR pre¬ 
sented in Ref. [16] is based^ on time-inversion symme¬ 
try. However, with the observation of charge-parity (CP) 
violation [17, 18], the charge-parity-time (CPT) theo¬ 
rem [19-21] suggests that time-reversal is no (exact) sym¬ 
metry either. Here, on the other hand, the argument 
for the parity SSR is constructed such that we rely only 
on Lorentz invariance via the spin-statistics connection, 
which happens to also be a central requirement for CPT 
symmetry. On the other hand, the SSR constraint im¬ 
posed in this way provides the essential ingredient to 
guarantee a meaningful quantum (information) theory 
based on anticommuting operators, as we will show next. 


V. SUPERSELECTION RULES & SYMMETRIC 
PURE STATE MARGINALS 

We shall now provide the main technical statement of 
this work, and its proof. 

Theorem. The marginals pa{b) = TrB(A)(||V'))(('0ll) 
of any bipartition {A\B) of a pure state ||'0)) in the 
fermionic Fock space have the same spectrum, if 
satisfies the SSR prohibiting superpositions of even and 
odd numbers of fermions, i.e., 

11^)) satifies SSR ^ spetr(pyi) = spetr(pB) V (A|i3). 

Proof. To show this, let us consider a pure state ll'i/'iten)) 
in an n-mode fermionic Fock space, where N = 
{1, 2,..., n} denotes the set of modes, and without loss 
of generality we have chosen H'i/'^en)) to be a superposi¬ 
tion of states with even numbers of fermions. The set N 
is then partitioned into the subsets M = {p,i \ pij S N : 
pi ^ Pj'Ai ^ j-,i, j = 1,2,... ,m < n} and = N\M, 
such that N = Mu and M C = 0. With respect 
to this bipartition, we may write the state || V'CTen)) ttie 


^ Indeed, Ref. [16] hints at the possibility of basing the proof on 
rotations instead of time-inversion, but this argument is not elab¬ 
orated on in [16] and was only later published in [22]. 
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pure state decomposition 

m 

llV'e^en)) = 70 || 0 )) + ^ 7 m. II 1 ^. )) II V'/lfodd )) ( 9 ) 

m 

+ E^AU..I|lM.))l|l.,))lK;..eve„)) + --- 

* J = 1 
j>i 

+ l|lMl))---l|lM™))IIV’E.M,..odd)) : 

where 7^., ■ • ■.^ <C, (ipevenHeven)) = 1. and 
without loss of generality we have here selected m to 

be odd. The states II V'E.A*fc.eve„(odd) )) contain only even 
(odd) numbers of excitations, and only in modes from 
the set . Any sign changes that may occur when 
rewriting a given state in such a decomposition can 
be absorbed into the y-coefficients. Adhering to the 
“outside-in” tracing rule of Eq. (5), we note that the 
state has been brought to a form, where the partial trace 
over is achieved by simply removing all projectors 
IIV’E.Mfc))((V'E.MJI pertaining to from the projec¬ 
tor on ||V’CTen))i without incurring any additional sign 
flips. On the other hand, if we trace over the modes in 
the set M instead, anticommuting the operators corre¬ 
sponding to modes in M towards the vacuum projector 
in the process, we may generate sign changes. However, 
for the superselected state, all the nonzero contributions 
to the partial trace over M are generated from elements 
such as 

TrM (1 7/., P 11 ^/^.odd)) ((V-ilf.odd 11 

= | 7 mJ^ lK.odd))(«.oddll • ( 10 ) 

There, the parity of the number of anticommutations to¬ 
wards II0)) from the left, is the same as the parity of the 
number of anticommutations towards ((0|| from the right. 
In other words, once the state has been brought to the 
form of Eq. (9), the partial trace can be carried out as if 
operating on a tensor product of Hilbert spaces, ie., as 
if Hn = Hm ®'Hm^ ■ la particular, this implies that the 
two reduced states on "Hm and which are isomor¬ 

phic to the corresponding m-mode and (n — TO)-mode re¬ 
duced fermionic states, respectively, have the same spec¬ 
trum. For any bipartition of the Fock space, a decom¬ 
position with this property can be found, although, in 
general, no decomposition exists that simultaneously ac¬ 
complishes the required task for all bipartitions at once. 
An example is presented in the Appendix. We hence con¬ 
clude that the SSR forbidding superpositions of even and 
odd numbers of fermions, guarantees that the spectra of 
the marginals for any bipartition are pairwise identical. 
An analogous argument applies if the initial state is a 
superposition of only odd numbers of fermions, or if m 
is even, and the proof therefore applies without restric¬ 
tion. □ 


VI. DISCUSSION 

Within quantum information theory and quantum 
computation, discussing problems in an abstract context 
has proven to be very useful. However, when attempt¬ 
ing a similar approach to a quantum information theory 
based on a fermionic Fock space one encounters difficul¬ 
ties. As we have shown, an unrestricted fermionic model 
features pure states with non-symmetric bipartitions. 
That is, pairs of reduced states across bipartitions need 
not have the same spectra, which is very problematic for 
the definition of entanglement. As we have shown, this 
problem is removed, when superpositions of even and 
odd numbers of fermions are forbidden. Nonetheless, 
the removal of an inconvenience appears to be a rather 
weak justification for the introduction of a restriction of 
generality within the abstract model. On the other hand, 
when placing the fermionic system within the physical 
framework of a relativistic quantum field theory, the 
SSR follows naturally from the requirements of Lorentz 
invariance and symmetry under rotations by 2tt. 

We hence argue that, in contrast to bosonic or 
qudit-based variants, any fermionic quantum informa¬ 
tion theory must be seen as (part of) a relativistic 
quantum field theory. Without such a physical context, 
a reasonable quantum information theory based on 
anticommuting operators can be obtained by adding the 
parity SSR as an axiom, but such an approach would 
seem to be rather ad hoc. This strong hint at the insep¬ 
arability of quantum information theory and the theory 
of relativity is rather surprising, but may provide deeper 
insight into constructions of quantum theory based on 
informational principles, see, e.g., Refs. [23-25]. More¬ 
over, our results provide fresh insight into the debate of 
entanglement in systems of indistinguishable particles 
(see, e.g., [35, 36]) in general, and questions of entangle¬ 
ment between fermionic modes [40] specifically. Finally, 
it will be of significant interest to see to what extent 
simulations of fermions, for instance in superconducting 
materials [41], or in graphene [42, 43], can capture the 
behaviour of superpositions of different particle numbers. 
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Appendix: Superselected fermionic pure state 
ineqnivalent to qubits state 

Here, we present an example for a pure state of four 
fermionic modes that satisfies the superselection rule 
(SSR) that forbids superpositions of even and odd num¬ 
bers of fermions. We explicitly compute the marginals 
of this state and show that the subsystem spectra match 
for any bipartition. We further prove that, nonetheless, 
this state and its marginals do not admit a consistent 
mapping to a four-qubit state (and its marginals). The 
most general pure state for four fermionic modes, labelled 
1, 2, 3,4, that only contains even numbers of excitations, 
may be written as 

IIV'^vD) = «o||0)) + ai 2 ||li))||l 2 )) + ai3||li))||l3)) 

+ ai4 II li))|| I 4 )) -I- Q!23 II l 2 ))|| I 3 )) 

+ 0^24 II l2))|| I 4 )) + a34 II l3))|| I 4 )) 

+ ai234||ll))||l2))||l3))||l4)), (A.l) 


be., operators pertaining to modes that are being traced 
over are anticommuted towards the projector on the 
vacuum state, before being removed. This prescrip¬ 
tion, which uniquely determines the marginals (see, e.g., 
Ref. [15]), is a direct consequence of the requirement that 
expectation values of local observables give the same re¬ 
sult when evaluated for the global state, or for the re¬ 
duced state, that is. 


■ (A.5) 

For the marginals of the bipartition (1, 2|3,4) we find 

P 1.2 = Tr3.4 (II ((4tL II) = PIT p!T + Pif Pif ’ 

(A.6a) 

P3.4 = Tri,2(||V'^^L))((^ivLll) =PlTP^^r+P3AP3A ^ 

(A.6b) 


where the reduced state density operators in the even 
and odd subspaces are given by 


where |aoP -f |ai 2 p -I- |ai3p + |ai4p + |a23p + |a24p + 
|a 34 p + |cn 234 p = 1- First, let us consider the biparti¬ 
tions into two subsets containing one and three modes, 
respectively. For instance, let us consider the bipartition 
(1|2,3,4). We may decompose the state HV’even)) into 
terms containing excitations in mode 1, and those which 
do not contain such terms, be., we write 

IIV'(eL)) = 71 ll<(>i))+7^1 110 ^ 1 )), (A.2) 

where Iqip -p Iq^^p = b\bi\\(l)i)) = ||(()i)), and 

bi llt/i^i)) = 0. Since ((^1 ||(('^i)) = 0, the marginals with 
respect to the bipartition (1|2,3,4) are then easily ob¬ 
tained as 

Pi= Tr2.3,4(||i/’ive„))((V’evLll) (A.3a) 

= l 7 in|li))((li|| + I 7 ^ini 0 ))(( 0 || , 

P2.3.4= Tri(||i/>W^))((V^(4Lll) (A.3b) 

= l7ini^i))((^i|| + l7^iPll<^^i))((<(>^i|| , 

where ||^i))((^i|| = Tri(||(/)i))(((/)i ||). Because HV'evL)) 
contains only even numbers of fermions, the same is 
true for ||^i)) and ||0^i)), whereas ||^i)) contains only 
odd numbers of fermions. Hence we may conclude 
that ((<^ 110 ^ 1 )) = 0, and it is thus easy to see that the 
marginals have the same spectrum. The same argument 
goes through for the bipartitions (2|1,3,4), (3|1,2,,4) 
and (4|1,2,3). 


PIT PIT = (l«oP + |a34n ||0))((0|| (A.7a) 

+ (|ai2p + |ai234n I|ll))||l2))((l2||((ll|| 

+ (aoCll2 + «34ai234) I|0))((l2 ||((li || + H.c. 

PifpT= (I«i3p + |ai4p) ||li))((li|| (A.7b) 

+ (|a23p + l«24p) ||l2))((l2|| 

(ai3a23 + «14a24) I|ll))((l2|| +H.C. 


for the subspace of modes 1 and 2. Similarly, for 3 and 4 
we obtain 


pi:rpi:r = (i«op + iio))((oii (A.sa) 

+ (|a34p + |ai234H I|l3))||l4))((l4||((l3|| 

+ (aoa34 + «i2a;i234) I|0))((l4 ||((l3 II + H.c. 

P3AP3f= (l«13r + l«23|") ||l3))((l3|| (A.Sb) 

+ (|ai4p + |a24p) ||l4))((l4|| 

(Q!l3a*4 -I- a23a24) II l3))((l4 II + H.c. 


For the superselected state the even and odd subspaces 
decouple. We may therefore compare the characteristic 
polynomials for the even and odd subspaces separately. 
A simple computation reveals that both p®™" PiT ^^id 
PsT PsT the characteristic polynomial 


Let us now turn to the bipartitions into pairs of modes. 
We use the tracing rule 

Trk{bl...bl bl\\0)){{0\\bkb,,...bT 

= II 0))((0 II V > (A.4) 


det{pTT pTT - ^1) = det(piyr pTT - Al) 

= — A(|aoP + |ai2p + WzT + |c4i234p) 

+ (aoai234)(cioai234)* + (ai2a34)(ai2a34)* 

— (aoai234)(ai2Q!34)* — (aoai234)*(cii2a34) ■ (A.9) 
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Similarly, for the odd subspace we find 

det{p°f p°f - Al) = det(p°;“ - Al) 

= A^ — A(|ai3p + |Q!i4p + |Q;23p + |Q^24p) 

+ (ai3a24)(Q!l3Q!24)* + (ai4a23) (c«l4Q!23)* 

— (ai3a24)(ar4a23)* — (ar3a24)*(ar4a23) ■ (A.10) 

We hence find that the eigenvalues of the marginals pi ,2 
and p 3^4 coincide. Some straightforward computation 
along the same lines conhrm that this is also the case 
for the bipartitions (1,3|2,4) and (1,4|2, 3). Now, the 
interesting aspect of this insight pertains to the fact that 
the four-mode system may not be consistently mapped 
to a four-qubit state. For the latter, the matching sub¬ 
system spectra would be guaranteed by the Schmidt de¬ 
composition. Recall that the consistency conditions (see 
Ref. [15]) for partial tracing demand that the numerical 
value of the expectation values of any “local” operator 
(in the sense of mode-subspaces) is independent of being 
evaluated for the overall state, or for the corresponding 
local reduced state. These consistency conditions then 
fix the relative signs of different contributions from ma¬ 
trix elements of the total state to the matrix elements of 
the reduced states. For the example state at hand, the 
resulting off-diagonal matrix elements of all two-mode vs. 
two-mode bipartitions are collected in Table I. 

Now, one wonders, whether the pure state ||V'evL)) 
and its marginals can be faithfully represented as a 
four-qubit state IV'evL) G (C^)®"‘. Here, we will call 
such a mapping faithful, if all the diagonal matrix ele¬ 
ments of I "06™!!) (0even | and its marginals with respect to 
the computational basis match the diagonal elements of 
||0’evL))((0evL|| with respect to the Fock basis. For the 
off-diagonal elements we impose a slightly weaker con¬ 
dition, i.e., that the absolute values of the off-diagonals 
(with respect to the respective bases) match. This corre¬ 
sponds to demanding that measurements in the Fock ba¬ 
sis are reproduced, and that the marginals have the same 
spectra. These conditions imply that a faithful mapping 
from the Fock basis of four fermionic modes to the com¬ 
putational basis of four qubits must be of the form 



II0)) 



0000) , 

(A.lla) 

ll))l 

12 )) 


g»012 

11100) , 

(A.llb) 

ll))l 

13 )) 


g»013 

11010) , 

(A.llc) 

ll))l 

14 )) 


g*014 

11001) , 

(A.lld) 

l2))| 

13 )) 


g*<i>23 

10110) , 

(A.lle) 

l2))| 

14 )) 


g*<i>24 

10101) , 

(A.llf) 

l3))l 

14 )) 


g*<i>34 

10011) , 

(A.llg) 

I 3 )) 

14 )) 


g4<i>i234 \ nn) . 

(A.llh) 



({ 0 ||p.,i||L))||li)) 

({ 1* II Pi,j II I7 )) 

1,2 

O0O12 + 03401234 

0130)3 -1- 0140)4 

3,4 

O0O34 -t 0120)234 

0130)4 -I- 0230)4 

1,3 

000)3 — 0240)234 

0140)4 — 0120)3 

2,4 

O0O24 — 0130)234 

0120)4 — 0230)4 

1,4 

000)4 + 0230)234 

—0120)4 — 0130)4 

2,3 

O0O23 -|- 0140)234 

0^120^13 + 0^24(^34 


TABLE I. The off-diagonal matrix elements of the reductions 
pij = Tr^i,j(||0ivL))((0even||) are shown. 

Performing the mapping of (A.11) for the state 

of Eq. (A.l), i.e., H^ivL)) i-T |0evL), and tak¬ 
ing the partial traces for the qubits as usual, 
we obtain the off-diagonal elements that are to 
be compared with those in Table 1. For instance, 
comparing | (0000 | Tr 3_4 (| ipitln ) ('0even I) 11100) I with 
I ((0|| Pi ,2 ||li))||l 2 )) I we get the condition 

= \a 0 al 2 + a34<ai234l • (A.12) 

Since this must hold independently of the values of 
Q!i 2 , Q! 34 , and 01234 , we arrive at 

012 + 034 — 00 ~ 01234 = 2 ni TT , (A.13) 

where ni G Z. Similarly, the other off-diagonals from 


Table I provide the conditions 

013 + 024 - 00 - 01234 = (2 712 + 1) TT , (A.14a) 

014 + 023 - 00 - 01234 = 2 773 TT , (A. 14b) 

023 + 014 - 013 - 024 = 2 774 77 , (A.14c) 

012 + 034 - 014 - 023 = (2 775 -h 1) TT , (A.14d) 

013 + 024 - 012 - 034 = 2 776 TT , (A.14e) 


with 772,3,4,5,6 G These conditions cannot all be met at 
the same time. This can be seen, for instance, by com¬ 
bining (A.13) with (A.14b), and comparing to (A.14d), 
which results in 

2(771 - 713 ) = (2 775 -I-1), (A.15) 

which cannot be satisfied, since the left-hand side is an 
even integer for all 771,773 G Z, while the right-hand side 
is an odd integer for all 775 G Z. We hence conclude that 
the state of Eq. (A.l) cannot be consistently mapped 
to a four-qubit state, even though it satishes the SSR, 
and despite the fact that for any of its bipartitions the 
respective marginals have the same spectra. 
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